Abstract. We report the results of a detailed analysis of 
the temperature structure of the X-ray emitting plasma 
halo of M 87, the cD galaxy of the Virgo Cluster. Using 
the MEKAL model, the data provide strong indications 
that the intracluster medium has a single phase structure 
locally, except the regions associated to the radio struc- 
tures. The deprojected spectrum at each radius is well 
fitted by a single temperature MEKAL model, except for 
the very central region (< 2 arcmin) which seems to be 
affected by the jet and radio lobe structure. The tempera- 
C^l , ture of the intracluster plasma is 1 keV at the center and 
j gradually increases to 2.5 keV at 80 kpc. We have also 
fitted spectra using the APEC code. Although the large 
changes of the strength of Ka lines causes a discrepancy 
^ , between the Fe-L and Fe-K lines for the APEC results, 
" the overall temperature structure has not changed. There 
ly-^ . is no sign of excess absorption in the spectral data. The 
single-phase nature of the intracluster medium is in con- 
flict with the standard cooling flow model which is based 
on a multi-phase temperature structure. In addition, the 
, signature of gas cooling below 0.8 keV to zero tempera- 
^-j- ' ture is not observed as expected for a cooling flow. The 
i gravitational mass profile derived from the temperature 
and density distribution of the intracluster gas shows two 
distinct contributions that can be assigned to the gravi- 
tational potential of the cD galaxy and the cluster. The 
, central temperature of the intracluster medium agrees well 
' with the potential depth and the velocity dispersion of the 
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cD galaxy. The latter result implies that the central region 



O ' of the intracluster medium is equivalent to a virialized in- 
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terstellar medium in M 87. 
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1. Introduction 

In the cores of many clusters of galaxies, X-ray imaging 
data show a highly peaked surface brightness profile (e.g. 
Fabian et al. 1981). The radiative cooling time in these 
regions is much less than a Hubble time. Without a heat- 
ing process, the gas cools to low temperature and results 
in a "cooling flow" (Fabian 1994 for a review). The mass 
flow rate, M, that is deduced in the standard cooling flow 
model, is approximately proportional to the radius. This 
implies that matter is deposited throughout the entire 
cooling flow region. It also implies that the gas in the 
cooling flow zone is "multi-phase" on scales small enough 
that the inhomogeneities have escaped the observation so 
far. 

ASCA and ROSAT observations confirmed the exis- 
tence of cooler gas in the cores of cooling flow clusters 
as expected (e.g. Allen & Fabian 1994, Ikebe et al. 1999, 
Ikebc 2001). The presence of an intrinsic absorber in ex- 
cess of the galactic absorption, as it might be expected 
to arise from the accumulation of cold gas in the cooling 
flows, was also inferred from the analysis of the ASCA 
spectra (e.g. Allen 2000; Allen et al. 2001). In fact, the 
fitting of the "multi-phase" cooling flow model to the ob- 
served spectra including temperature phases that cool be- 
low the X-ray emitting temperature regime, requires the 
inclusion of excess absorption to satisfactorily reproduce 
the observed spectra. It has also been argued, that inclu- 
sion of excess absorption leads to values of the mass flow 
rate as obtained from imaging data, M\ in good agreement 
with those obtained from spectral analysis, Ms (Allen 
2000). 

However, without assuming excess absorption, the val- 
ues of Ms are systematically lower than M\ (e.g. Ikebe 
et al. 1999; Makishima et al. 2001). Instead of the cool- 
ing flow model, a two temperature spectral model can 
well fit ASCA spectra of the cooling flow clusters. The 
temperatures of the softer component are about half of 
the hotter component (Ikebe 2001). The central tempera- 
tures of some cD galaxies obtained with the ROSAT PSPC 
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are close to those obtained from normal elliptical galax- 
ies with the same stellar velocity dispersion (Matsushita 
2001). Since cooling flow clusters always possess central 
dominant galaxies, the cool components may reflect the 
potential structure, rather than cooling gas (Ikebe et al. 
1999, Makishima et al. 2001, Ikebe 2001). 

Recently, it was discovered with the RGS instrument 
onboard XMM-Newton that there is little X-ray emission 
from a component with a temperature below a certain 
lower cutoff value that differs from object to object (e.g. 
Tamura et al. 2001; Kaastra et al. 2001). However, with 
the RGS one observes only the very central part of the 
cluster. 

M 87 is the cD galaxy of the nearest rich cluster of 
galaxies, the Virgo Cluster. It is very luminous in X-rays, 
and is suggested to have a "cooling flow" with M of about 
10 Moyr" 1 (Stewart et al. 1984; Fabian et al. 1984). M 87 
also hosts a central active galactic nuclei (AGN). The ra- 
dio emission is complex and there are two strong lobe 
structures (e.g. Bohringer et al. 1995). An enhancement of 
the X-ray emission around the lobes was discovered with 
the EINSTEIN HRI (Feigelson et al. 1987) and ROSAT 
PSPC (Bohringer et al. 1995) 

In addition to the XMM RGS results mentioned above, 
the XMM-Newton observatory offers the possibility to 
perform a very detailed spectral and spatial analysis of 
the cooling flow area. This is in particular important 
since the RGS spectrum taken at the very central part of 
the intracluster medium (ICM) may be strongly affected 
by resonance line scattering and thus the interpretation 
of the results requires the combined analysis of spectra 
taken at different radii from the center. Early results of 
the XMM observation of M 87 have already been pub- 
lished (Bohringer et al. 2001; Belsole et al. 2001). The pro- 
jected spectrum is fitted with a single temperature model 
(Bohringer et al. 2001). The abundance drop at the central 
region is interpreted as resonance scattering. The detailed 
spectral analysis on the jets and the radio lobe structure 
is shown in Belsole et al. (2001). 

In this paper, we report on the more detailed temper- 
ature structure using a deprojecting analysis, and apply a 
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Fig. 1. EPN image of M87 in the energy range of 0.2 to 
10.0 keV. The effect of the out-of-time events is corrected. 
The overlayed contours correspond to the image divided 
by a azimuthally symmetric with the same radial profile, 
smoothed with 20" Gaussian. The contour levels are 0.8, 
0.9 (white), 1.1, 1.2, 1.3, 1.4, and 1.5 (red). The regions 
surrounded by the blue contours are excluded from our 
spectral analysis as described in §3.3. 

cooling flow model to the spectral analysis. We employ a 
distance of M 87 of 17 Mpc. We adopt for the solar iron 
abundance the 'photospheric' value, Fe/H — 4.68 x 10 -5 
by number (Anders & Grevesse 1989). Unless otherwise 
specified, we use 90% confidence error regions. 

2. Observation and data preparation 

M87 was observed with XMM-Newton on June 19th, 2000. 
The thin filter is used for both detectors. 

2.1. EPN 

Figure 1 shows the X-ray image of M87 obtained with 
the EPN. The contribution from the out-of-time events is 
subtracted following the standard analysis package (SAS). 
The image divided by an azimuthally symmetric image 
with the same radial profile is overlayed in Figure 1. The 
X-ray halo is very smooth and almost spherically symmet- 
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ric within 20%, except for the two inner radio lobe regions 
and the first knot in the jet. A luminous point source near 
the south edge of the detector is excised in the further 
analysis. 

We have selected single pattern events. The back- 
ground spectrum was calculated for all spectra by inte- 
grating the Lockman-Hole data (observed in the Revolu- 
tion 70,71,73,81 for the EPN) in the same detector re- 
gions. In the Lockman-Holc field, sources with a luminos- 
ity higher than that of the excluded sources from the M 87 
region are excised. The flare events in the background data 
were screened in the following way. We made two count 
rate histograms of bin width of 200s, for the corner region 
in the whole energy band, and for the whole detector for 
above 10 keV. Then, we fitted each histogram with a gaus- 
sian, and selected the time within 3 a of the mean for both 
histograms. The total exposure of the background data is 
95ks. We have also screened the M87 data in the same 
way. The exposure times of the EPN are 30ks, 28ks, 27ks, 
30ks for the Chips 1-3, 4-6, 7-9, and 10-12, respectively. 

The spectra of each CCD chip are accumulated within 
a ring, centered on M 87. Then, we have subtracted the 
background and corrected the effect of the vignetting and 
exposure, and summed the spectra within annular regions. 
The contribution of out-of-time events is subtracted using 
SAS. 

2.2. EMOS 

We have selected events with pattern between to 12. The 
background spectra are accumulated from the Lockman- 
Hole data observed in Revolution 70. The background flare 
events are screened in the same way for the EPN. The 
exposure times of the background and data are 33 ks and 
40 ks, respectively. Out-of-time events are not subtracted 
since the contribution is negligible. 

2.3. Spectral Deprojection 

Deprojected spectra are calculated by subtracting the con- 
tribution from the outer shell regions for all spectral com- 
ponents, assuming the ICM is spherically symmetric as 
done by Nulsen and Bohringer (1995). Within a each shell, 
the spectrum per unit volume is assumed to be the same. 
In order to calculate the deprojected spectrum of the out- 
ermost (N th ) shell, we subtracted the emission of the out- 
side of the N th shell using the brightness profile derived 
from the ROSAT All Sky Survey (Bohringer et al. 1994) 
assuming that spectral shape outside the N th shell is the 
same with that of the N th shell. The deprojected spec- 
trum of the i th shell is then calculated by subtracting the 
contributions from i + 1 th to N th shell and from outside 
of the N th shell from the annular spectrum of the corre- 
sponding radius. 
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Fig. 2. Annular spectra of the EPN (black) and EMOS 
(EMOS1+EMOS2; gray) for r=5.6-8', fitted with a same 
single temperature MEKAL model. Bottom panels show 
data-to-model ratios. 

3. Spatial distribution of temperature 
components 

3.1. Radial temperature profile 

We employ epn_fs20_sY9_thin.rmf from March 2001, 
ml_thinlv9ql9t5r5_all_15.rsp from June 2001 for the EPN 
and the EMOS response files, respectively. The spectral 
analysis uses the XSPEC.vl 1.0.1 package. 

Figure 2 shows annular spectra for the radial range 
r=5.6-8' of the EPN and the EMOS, fitted with a same 
single temperature MEKAL model (Mewe & Gronenschild 
1985, Mewe et al. 1986, Kaastra 1992, Liedahl et al. 1995) 
but different normalization. Here, r is the the projected ra- 
dius. The EMOS gives about 10% higher normalizations 
than the EPN. The response matrices of the EPN and 
EMOS arc mostly consistent. Below 0.5 keV, EPN data 
is lower than the EMOS data, and there arc small dis- 
crepancies around the tail of the Fe-L lines at 0.8-1.0 keV 
and around the instrumental Si edge at 2 keV. Since the 
EMOS gives a hydrogen column density, Ah, consistent 
with that of ROSAT (c.f. §6), we fitted to the EPN and 
EMOS data in the spectral range of 0.5 to 10 keV and 
0.3 to 10 keV, respectively. For R > 8', the energy range 
between 7.5 and 8.5 keV of the EPN spectra is ignored, 
because of strong instrumental emission lines. 

We have fitted the spectra for the region with radius 
R > 0.5' with a MEKAL model with photoelectric absorp- 
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Fig. 3. Dcprojected radial profile of the temperature for 
the EPN (gray crosses) and the EMOS (gray diamonds), 
obtained using the single temperature MEKAL model. We 
also plotted the same result but excluding the jet and ra- 
dio lobe regions (Figurel; §3.3) for the EPN (black crosses) 
and the EMOS (black diamonds). The solid line and the 
dashed line corresponds to the best-fit power-law relation 
(excluding the lobe regions) for the EMOS and EPN, re- 
spectively. 

tion. Here, R is the three dimensional radius. For the spec- 
tra within 0.5', a power-law component, with the same ab- 
sorption, Ah, is added. When fitting the spectra between 
0.12' and 0.5', we fixed the index of the power-law com- 
ponent to the best-fit value obtained from the spectrum 
within 0.12' and normalized it using the Point Spread 
Function. We grouped several elements and constrained 
the elements in each group to have a common abundance. 
The first group is C, N, and O, and the second group is Fe 
and Ni. The other elements are determined separately. We 
also fitted the spectra with a sum of two MEKAL mod- 
els, where the two components were constrained to have 
common abundances. 

We have summarized the results in Tabic 1. 

Figures 3 and 4 show the radial profile of the temper- 
ature obtained by the single temperature model, and the 
two temperature model, respectively. The temperatures 
derived from the EPN and EMOS are mostly consistent. 

Within R < 6' and R < 13.5' of the EMOS and EPN 
respectively, the two temperature model gives significantly 
lower \ 2 values than the single temperature model. The 
temperature of the cooler component is almost constant 
at ~ 1 keV, which is the same as the central temperature 
obtained by the single temperature model fit. That of the 
hotter component show a nearly constant value at R < 2', 
whereas for R > 2' it increases gradually. In Figure 5, 
we have plotted the ratio of the emission measure of the 
cool to the hot component, assuming the temperature of 
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b: Ratio of emission measure of the cooler and hotter component, when the temperature of the former is fixed to 1 keV. 
Table 1. Result of spectrum fitting of the deprojected spectra (whole region) 



the former to be 1.0 keV. The fraction of the cool com- 
ponent is 20% at R=0.5' and drops below 1 % beyond 6' 
for the EMOS. Beyond 4', the EPN tends to have 1-2% 
larger values of the contribution of the cooler component 
than the EMOS. This difference reflects the discrepancy 
at 0.8-1 keV between the EPN and the EMOS shown in 
Figure 2. This energy band corresponds to an instrumen- 
tal low-energy tail of the Fe-L lines and the effect for the 
EPN is stronger than that of the EMOS. The EPN has 
an additional uncertainty due to the position dependent 
energy resolution. Considering that the EMOS spectra of 
R > 6' at 0.7-1.1 keV are well fitted by the single tem- 
perature MEKAL model fit, the larger amount of 1 keV 
component observed by EPN may be caused by the uncer- 
tainty of the response matrix of the EPN. Thus, there is 
a systematic uncertainty of 1-2% in the fraction of the 1 
keV component. In summary, within R < 6', we detected 
at least two temperature components and outside R > 6', 
the upper limit of the second component is less than a few 
% even using the EPN. 



ASCA spectra of M 87 are also well fitted with a sum 
of two Raymond-Smith models (Raymond & Smith 1977). 
The temperature of the two components are 3.0 and 1.3 
keV, and the ratio of the emission measure range from 0.2 
to 2 within 10' (Matsumoto et al. 1996), although the re- 
sults of MEKAL model fit with ASCA spectra are mostly 
consistent when considering the difference in spatial res- 
olution (Finoguenov & Jones 2000; Shibata 2000). The 
contribution of the cooler component plotted in Figure 5 
is much smaller than that derived from R-S fits. 

3.2. Azimuthal temperature structure 

In a second step, in order to study the spatial distribu- 
tion of the temperature components in more detail, the 
EMOS spectra are accumulated within sector regions as 
shown in Figure ^|. The spectral fits used two temperature 
MEKAL model. Although the spectra fitted in this sec- 
tion are not the deprojected ones, we can constrain the 
spatial distribution of the temperature structure. Figure 
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Fig. 4. The radial temperature profile obtained by fitting 
the spectra with a two component MEKAL model. The 
meanings of symbols is the same as in Figure 3. The solid 
(EMOS) and dashed line (EPN) corresponds to the best- 
fit power-law ration of the single component MEKAL fit, 
excluding the radio lobe regions. 



| shows EMOS images of the Fe-L (0.6-1.2 keV) and hard 
(2-10 keV) energy band with the temperatures obtained 
for the two components. It also shows hardness ratio maps 
for the energy band ratios 0.6-1.2keV/2-10keV and 1.0- 
1.2keV/0.8-1.0keV. Although the latter is the most sensi- 
tive to the 1 keV components (Figure 7), the statistics is 
not sufficient in the outer regions. The ratios of the emis- 
sion measure of the cool to hot components are overlayed 
on the hardness ratio maps. 

The hard band image is nearly spherically symmetric 
and no enhancements of the X-ray emission associated to 
the radio structure is seen. The temperature of the hot 
components is also almost spherically symmetric within 
10%. It tends to be slightly lower in the region around the 
radio lobes. 

In contrast, the spatial distribution of the cool com- 
ponent is far from symmetry, although its temperature is 
determined to have an almost constant value, ~ 1 keV. 
Belsole et al. (2001) found that the softer regions agree 
well with the radio structures. Within R <1\ the 1 keV 
component is detected from all the spectra, but its con- 
tribution ranges from several % to 30 %. In the radial 
zone R=l-2', where 4 sectors are located outside the ra- 
dio lobes, we find a signature of the 1 keV component in 
two of these sectors but it is absent in the other two. Out- 
side 2', we have only detected the 1 keV component from 
the two sectors containing the radio lobes. 




R (arcmin) 

Fig. 5. The ratio of emission measure of the cool and hot 
temperature component, assuming the temperature of the 
former to be 1.0 keV. The meanings of symbols is the same 
as in Figure 3. The discrepancies between the EPN and 
EMOS may be caused by instrumental problems in the 
EPN (see text). 

3.3. Excluding the radio lobe regions 

For the further analysis, we have filtered regions with ex- 
cess emission related to the jet and radio lobes. We used a 
spatial filter, excluding those regions where the brightness 
is larger by 15 % than the azimuthally averaged value (Fig- 
ure 1, ^). The softer regions are not fully excluded within 
2', since it is difficult to filter whole the 1 keV component 
within 2' due to its complicated structure and the limited 
spatial resolution of the XMM telescope. 

Assuming spherical symmetry, deprojected spectra are 
calculated. We also made deprojected spectra with larger 
shell width than those in Table 1. We accumulated spectra 
within annular regions again and subtracted the contribu- 
tion from the outer shells with thin shells width obtained 
before. We then fitted the spectra with the single temper- 
ature and two temperature MEKAL model, in the same 
way as in §3.1. In the following sections, we will only use 
the spectra filtered the regions with enhance emissions. 

The results are summarized in Table 2. After filtering 
the regions with excess emissions, the emission measure 
of the 1 keV component is reduced (Figure 5), and the 
single temperature model gives reasonable fits at R > 2' 
(Table 2; Figure [?]). The remaining 1 keV components of 
the EPN beyond 2' may be due to the uncertainty of the 
instrumental low-energy tail of the Fe-L lines. The radial 
temperature profile from the single temperature MEKAL 
model fit is approximated as kT = 1.52r^^i?°- 23 -n.ol keV 
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Fig. 6. Hard (2-10 keV; left panels) and Fe-L blend (0.6-1.2 keV; middle panels) EMOS images, overlayed by the 
hard and soft components temperatures of the two MEKAL component model fit to the projected sector regions, 
respectively. The 90% uncertainty of temperatures of the cool and hot components are typically 0.1 and 0.03 keV, 
respectively. The top panels and bottom panels are 16'xl6' and 4'x4' regions, respectively. The right panels show 
the hardness ratio maps of 2-10keV/0.6-1.2keV (top) and 1.0-1.2keV/0.8-1.0keV (bottom), overlayed by the ratio of 
the emission measure of the soft to hard components in unit of percent. White contours corresponds to the regions 
excluded from spectral analysis after §3.3. 



for the EMOS and kT = 1.51±j];^i? a21 -o "i keV for the 
EPN (Figure 3), where R is in units of arcmin. 

Within R < 2', the two temperature model still gives 
better fits than the single temperature model, because 
of the complex structure of the 1 keV within the re- 
gion. When we adopted the hotter component tempera- 
ture within 2', the radial temperature profile is approxi- 
mated as kT = 1.69±g;gf(l + (R/1.&)) - 115 -™™ keV and 

kT = 1.68±g;gg(l + (iJ/1.6'))°' 100 -» : °" keV for the EMOS 
and EPN, respectively. 

There are some residuals at ~1.2 keV, which are com- 
mon to all the deprojected shells. Any multi-temperature 
model cannot explain the residual, because the continuum 
and the Fe-K line are fitted well with a single MEKAL 
model. When the Ni abundance is assumed to be a factor 
of 2 larger than Fe in solar units, the residuals become 
smaller (Finoguenov et al. 2001; Matsushita et al. 2002). 
This structure may also be due to a problem the Fe-L 
atomic data, since the APEC code (Smith et al. 2001) do 
not need such a high Ni abundance. 



These results indicate that most of the space is oc- 
cupied by a component which is nearly spherically sym- 
metric. Within relatively small regions associated to the 
radio structures, there is an additional component whose 
temperature is about 1 keV. In the following sections, we 
will study the temperature structure in detail using the 
filtered spectra. 

4. The abundance profile of Fe 

In spectral fitting the derived temperature and element 
abundances are partly dependents on each other. This is 
particularly important for the Fe abundance, where the 
spectral features of K- and L-shcll lines with their tem- 
perature dependent relative importance are fitted simul- 
taneously. A more detailed report on the abundances is 
given in Paper II (Matsushita et al. 2002) and Finoguenov 
et al. (2001), however we briefly discuss here the results 
on the Fe abundances in relation to the derived tem- 
perature structure. As derived from the projected spec- 
tra (Bohringer et al. 2001; Finoguenov et al. 2001), the 
Fe abundance shows a strong negative gradient (Figure 
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b: Ratio of emission measure of the cooler and hotter component, when the temperature of the former is fixed to 1 keV. 
Table 2. Result of spectrum fitting of the deprojected spectra (regions with the enhanced X-ray emission associated 
to the radio structures are excluded) 



H ). The two temperature MEKAL fit gives significantly 
larger abundances than the single temperature MEKAL 
fit, although the contributions of the cooler component are 
small. For example, at i?=0.5-l', including of 10 % of the 
cooler component in units of emission measure changes Fe 
abundances by a factor 2. The reason for this is the change 
of the temperature of the hot component, when the cold 
component is added, e.g. in the ring taken as an example 
the temperature of the main component increases from 
1.47 to 1.72 keV. For a fixed Fe line feature this leads to a 
higher abundance required in the fit. As a result, the cen- 
tral abundance drop which is seen in the single MEKAL 
fit is not seen in the result of the two MEKAL fit. This will 
be discussed further in paper II (Matsushita et al. 2002). 



5. Detailed temperature structure 

5.1. Temperatures obtained from the Fe-L and the hard 
energy band 

There are some problems in the determination of the tem- 
perature through X-ray spectral fitting. The first problem 
is resonant line scattering, which should be important in 
the core of M 87 (Shigeyama 1998, Bohringer et al. 2001). 
The second problem is the uncertainties in the Fe-L atomic 
data (e.g. Masai 1977; Matsushita et al. 2000). The major 
uncertainty is the ionization and recombination rates of 
Fe, which may change resulting temperature by 20-30%. 
(Masai 1997; Matsushita et al. 2000). In order to avoid 
these problems, we exclude the Fe-L region below 1.6 keV 
and fit the spectra again. In addition, to compare the tem- 
perature obtained to those from the Fe-L lines, we also fit- 
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Fig. 7. Deprojected spectra of the EPN (lobes excluded). 
The spectrum at R < 0.12' is fitted with the double com- 
ponent model, consisting a MEKAL model and a power- 
law component and those of R=0. 5-1.0' and 1.0-2.0' are 
fitted with a two component MEKAL model. Each com- 
ponent is plotted in dashed lines. The other spectra are 
fitted with a MEKAL model. The bottom panels show 
residuals of the fit. 



ted the data between 0.7 to 1.3 keV with a single temper- 
ature MEKAL or Raymond-Smith model (hereafter R-S 
model). The R-S model is used for the comparison with 
the ASCA results. 

The result is summarized in Table 3 and Figures ^|. 
The derived temperatures using the hard energy band, 
fcTkard, of the EPN and EMOS are consistent with each 
other (Table 3). At R > 2', they are consistent to those 
obtained using the whole energy band and within i?=0.35- 
2' they are consistent to the temperatures of the hotter 
component obtained in §3.1 and 3.3. This indicates that 
within i?=0.35-2', the 1.7 keV component really exists 
and is the dominant phase in the ICM. 

The results for the Fe-L region using the MEKAL 
model, /cTmekalo are consistent or slightly higher than 
fc^hard at J? > 1' because of the residual structure around 
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Fig. 8. Deprojected Fe profile of the EPN (crosses) 
and EMOS (diamonds) fitted with a single temperature 
MEKAL model (black) and the two temperature MEKAL 
model (gray). 



~ 1.2 keV. Within R < 1', /cT M ekal is lower than fcT har d 
due to the contributions from the 1 keV components. 

The R-S model gives significantly lower temperatures 
than fcThard by ^20%. As a result, when fitted with a R-S 
model, there remain residuals in the hard energy band. To 
explain the residuals, we need an additional R-S compo- 
nent with a high temperature, ~ 3 keV. This is the reason 
of the comparable emission measure of the two tempera- 
ture R-S components obtained in the ASCA results (Mat- 
sumoto et al. 1996) 

5.2. Temperature obtained from line ratios 

The ratio of He-like K a to H-like K a line strengths is 
sensitive to the gas temperature. Since the energy resolu- 
tion of the EMOS is better than that of the EPN, we used 
only EMOS data here. We fitted the spectra within the en- 
ergy range of 1.8-2.1 keV, and 2.3-2.7 keV with a MEKAL 
model. We denote the temperature obtained from the for- 
mer energy range as fcTs; and from the latter as fcTg . The 
contribution of Fe-L lines at the energy of the He-like K a 
line of Si is considerable, and the He-like line contains sev- 
eral strong components. The result is summarized in Table 
4. 

Figure ^ also compares the temperatures, fcTsi and 
fcTg. These temperatures agree well within 10% with the 
deprojected temperature of the corresponding shell, al- 
though kTsi may have some uncertainty, since around the 
Si He-K lines, Fe-L and Ni-L lines are also important. Fur- 
thermore, there is an instrumental Si edge in the response 
matrix and a strong instrumental Si line, which may cause 
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a: kT obtained by fitting the spectra above 1.6 keV. 

b: kT obtained by fitting the spectra between 0.7 to 1.3 keV. 

Table 3. Result of spectrum fitting of the deprojected 
spectra (radio lobes excluded) 
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a: degrees of freedom 

Table 4. Result of spectral fit of EMOS data, using the 
energy band around the Si and S lines with a MEKAL 
model 
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Fig. 10. Deprojected EMOS spectrum within R=l-2.8' 
in the energy range of Si and S lines, fitted with a MEKAL 
model (black line) or 1.4 keV + 3.3 keV R-S model (gray 
line). 
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Fig. 9. Deprojected radial profile of the temperature ob- 
tained with EMOS, by fitting the spectra of 0.7-1.3 keV 
(black diamonds), above 1.6 keV (black crosses), 1.8-2.1 
keV (gray crosses) and 2.3-2.7 keV (gray diamonds) with 
a MEKAL model. The results obtained by R-S model us- 
ing 0.7-1.3 keV are also shown (dotted diamonds). The 
solid line corresponds to the best fit regression line for 
the MEKAL model using the whole energy band of the 
EMOS. 



the discrepancy between the EPN and EMOS at 2 keV 
shown in Figure 2. 

We also compare the model and data spectra of the 
Si and S region in Figure [10. The two temperature R- 



S model with the same abundance ratio, cannot fit the 
Si and S line ratios, while a single temperature MEKAL 
model fits the lines well. In order to explain the observed 
line ratios with the two component R-S model, the harder 
spectral component should have an abundance of Si and 
S a factor of 3 ~ 5 larger than the softer component. 

Since the deprojected data have larger uncertainties, 
we have also used the S line ratios of the annular spectra. 
The strength of the Ka lines are obtained by fitting the 
annular spectra within 2.2-2.8 keV with a bremsstrahlung 
and two gaussians. We calculated the projected S line ratio 
using the deprojected temperature and density (§7) pro- 
files and compared with the observed ones. We used two 
temperature relations derived in §3.3. One is the best fit 
relation of the single MEKAL model fit and the other re- 
lation is derived when we adopted the hotter component 
temperature of the two temperature MEKAL fit within 
R < 2'. Outside 1', the observed ratios agree well with the 
calculated profiles from single temperature model fit (Fig- 
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Fig. 11. Projected line ratios (EMOS) between Ka of 
hydrogen and helium like S (closed circles). Errors cor- 
respond to 68% confidence level. The solid and dashed 
lines represent the results using the observed temperature 
profile from the single temperature MEKAL model and 
that of the hotter component from the two temperature 
MEKAL model, respectively. The assumed radial temper- 
ature profiles are written in the figure. We also plotted the 
profile for temperatures shifted by 10% (dotted line) 



ure [Tl]) . Comparing these results with a profile calculated 
using a temperature profile shifted by 10% to the data, we 
can conclude that the temperature profile inferred from 
the observed S line ratios is consistent within several % 
with the deprojected temperature profile. Within 1', the 
observed ratios are close to with those calculated using 
the hotter component temperature of the two tempera- 
ture MEKAL fit, but smaller than those from the single 
temperature fit. This result indicates that S lines are dom- 
inated by those of the hotter component within 2'. 

We then compared the strength of the Fe-L and Fe-K 
lines as follows. As for the S line ratios, the line brightness 
profile of Fe-K is calculated using the deprojected results. 
Within 2', we adopted the Fe abundance obtained from 
the two temperature model fit. The Fe abundance from the 
deprojected data is mainly determined from the Fe-L lines. 
For a given temperature, this analysis gives a comparison 
between the Fe-L and Fe-K spectral signatures. Outside 
0.5', the observed profile of Fe-K is only 20% larger than 
the calculated profile (Figure |l^). Since the Fe-L atomic 
data may have uncertainties of 20-30% (e.g. Masai et al. 
1998; Matsushita et al. 2001), the difference may be due to 
the uncertainties in the Fe-L atomic data. If the difference 
is caused by the temperature structure, it is explained by 
only 10% larger temperatures. 
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Fig. 12. The projected radial profile of the line bright- 
ness of Fe-Ka of the EMOS (open circles) and the EPN 
(crosses). The solid and dashed lines represent the results 
using the temperature profile from the single temperature 
MEKAL model and that of the hotter component from the 
two temperature MEKAL model, respectively. The dotted 
lines corresponds to 10% higher temperatures. Errors rep- 
resent the 68 % confidence level. 



In summary, excluding the regions associated to the 
radio structures, at R>1', fcThardi kTgi and fcTg of a given 
deprojected spectrum agree within 10%. /cTmekal also 
agrees well with kT^^, although fcTp e -L depends on the 
plasma emission code (e.g Masai 1997). These results in- 
dicate that at least at R > 1', the ICM at a given radius 
is dominated by a single temperature component whose 
temperature is 1.7 keV at R—V and 2.5 keV at i?=15'. 
The consistency between Fe-L and Fe-K also supports the 
single phase temperature structure. The difference of these 
temperatures within 1' reflects that the contribution of the 
1 keV component become significant. 

Moreover, the temperature obtained from the MEKAL 
model seems to be reasonable at least in the obtained tem- 
perature range, when ignoring the residual structure at 1.2 
keV. 



5.3. Spectral fits with the APEC plasma code 

A new plasma code, APEC (Smith et al. 2001), is now 
available in XSPEC. The residual structures around the 
Fe-L emission that occurred in the fits with the MEKAL 
code disappeared with the use of the APEC code for Cen- 
taurus clusters (Sanders & Fabian 2001). 

In this paper, we used the version 1.0 of the APEC in 



the version 11.0.1 of the XSPEC. In figure 13, we show 



ratios of a single temperature plasma with solar abun- 
dance calculated with the APEC and MEKAL models, 
when observed with the EMOS. Within the energy reso- 
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Fig. 13. Ratio of theoretical spectra of MEKAL to APEC 
code for a plasma temperature of 1 keV (black solid line), 2 
keV (gray solid line) and 4 keV (black dashed line) when 
observed with the EMOS. Abundances are assumed to 
have solar values. 



lution of the CCD detectors, the dominant change is not 
in the Fe-L lines, but in the Ka lines. Using the APEC 
code, the strength of the Ka lines of hydrogen like ions, 
which arc the most fundamental lines, decrease by several 
tens of %, and the temperature dependence also changes. 
The 6.7 keV Fe-K lines decrease, by a factor of 2 at 1 keV 
and 30% at 4 keV. The changes of the Ka lines of helium 
like ions of Si, S, Ar and Ca are smaller. 



Figure 14 shows a representative EMOS spectrum of M 



87 fitted with the MEKAL and APEC model. As already 
shown, the single temperature MEKAL model can well 
fit the spectra including the strength of Fe-L, Fe-K and 
S lines. The APEC model gives a better fit for the Fe- 
L/Mg-K structure around 1.2-1.5 keV and the Fe-L/Si-K 
structure around 1.8 keV. However, the single temperature 
APEC model cannot fit the Fe-K and Fe-K, as well as the 
ratio of Ka lines of S simultaneously, although the fit to 
the Fe-L complex and the continuum is good. 

Using the APEC code, we fitted the deprojected spec- 
tra of the EMOS in the same way as described in §3.3, 
and §5.1-2. Figure [l5| and Tables 5 summarize the tem- 
peratures derived with the APEC code. Due to the incon- 
sistency of the S line ratios and the strength of the Fe-K 
line, the reduced x 2 is systematically worse at R > 2' com- 
pared to those obtained from the MEKAL model fit. The 
temperatures obtained from the whole energy band, the 
Fe-L complex and the hard energy band using the APEC 
code are consistent with those obtained using the MEKAL 
model. The derived Fe abundances, which are mainly de- 
termined by Fe-L lines, are also consistent with those from 
MEKAL model (figure [l7|). In contrast, the temperatures 
from Si and S energy band are 10% larger, due to the 
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Fig. 14. The deprojected spectrum of EMOS within 
#=4-8' (crosses) fitted with a MEKAL model (thin black 
lines), an APEC model (bold black lines) and a two tem- 
perature APEC model (bold gray lines). The bottom pan- 
els shows the contributions to x 2 . 



changes of the strength of the Ka lines of hydrogen like 
ions. 

The projected radial profile of the S line ratios are also 
compared with the calculated profile from the deprojected 
data (Figure [l6]). Although the deprojected temperatures 
cannot reproduce the observed line ratios, the difference 
is only 10% in temperature. Therefore, as in the case with 
MEKAL model, temperatures derived from whole energy 
band, the Fe-L, the hard band Si, and S energy band agree 
within 10%. 

However, there is a discrepancy between the Fe-L and 
Fe-K strength. Figure 17 also shows the Fe abundances 
obtained from the deprojected spectra above 2 keV, fitted 
with a single temperature MEKAL and APEC model. We 
denote the derived Fe abundances from > 2 keV, as Fe- 
Kmekal and Fs-KapeCj and those from the whole energy 
band as Fe — Lmekal and Ps-Lapec for the MEKAL and 
APEC codes, respectively. Outside R > 1', Fe-KMEKAL , 
Fe-LMEKAL, and Fe-LAPEC agree well each other. Within 
R < 1', Fe-KMEKAL are consistent with Fe abundances 
obtained from the two temperature MEKAL model. How- 
ever, Fe-KMEKAL are systematically larger by ^60%. The 
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Table 5. Result of EMOS spectrum fitting of the deprojected spectra using a APEC model (lobe excluded) 
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Fig. 15. Deprojected radial profile of the temperature ob- 
tained with EMOS, by fitting the spectra of whole energy 
band (bold black crosses), 0.7-1.3 keV (black diamonds), 
above 1.6 keV (black crosses), 1.8-2.1 keV (gray crosses) 
and 2.3-2.7 keV (gray diamonds) with a APEC model. 
The solid line corresponds to the best fit regression line 
for the MEKAL model using the whole energy band of the 
EMOS. 
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Fig. 16. Projected line ratios (EMOS) between Ka of hy- 
drogen and helium like S (closed circles). Errors corre- 
spond to 68% confidence level. Solid line and dashed lines 
represent the results from deprojected temperature pro- 
files using APEC model, fcT w h o i and /cTh ar d, respectively. 
We also plotted the profile for temperatures shifted by 
10% (dotted line) 



strength of the projected Fe-K lines are almost a factor of 
2 larger than the values calculated from the deprojected 
results (Figure |lj). In order to obtain the observed Fe- 



K line profile, the temperature should be shifted by 30% 
compared to the deprojected temperatures. 
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Fig. 17. Deprojected profile of Fe abundances obtained 
from the whole energy band (diamonds) and above 2 keV 
(crosses), using the single temperature MEKAL model 
(black), APEC model (gray) and the two temperature 
MEKAL model (dotted). 

Adding another temperature component, the discrep- 
ancy of Fe-K lines become smaller (Figure |l4j). However, 
the x 2 is still worse than the single temperature MEKAL 
model, and the radio of S lines cannot be well fitted. 

In summary, the changes in the version 1.0 of the 
APEC code for the Lya lines (H-like ions) are giving a 
much less consistent picture than the results from the 
MEKAL code. Therefore, we are not sure that the new 
atomic data included in APEC constitute an improve- 
ment. 



6. Intrinsic absorption 

The hydrogen column density, TVh, derived from the single 
MEKAL fit is almost constant at R > 2' (Table 2; Fig- 
ure |is| ). As summarized in table 2 and 5, the ATjj derived 
from the single APEC model is consistent with those from 
the single MEKAL model. The average value at R > 2' 
is 1.5 (1.2-1.7) xl0 20 cm- 2 and 1.3 (1.0-1.5) xl0 20 c m - 2 
for the EPN and EMOS, respectively. Within 2 arcmin, 
the two phase MEKAL model fit gives a consistent value 
with those at R > 2', although the single MEKAL model 
fit gives a higher A^h value by several xl0 20 cm~ 2 . The 
response matrices of the EPN of September 2000 and 
November 2000 give systematically different values by 1 
~ 2 xl0 20 cm~ 2 . Therefore, there may be a systematic 
uncertainty of a few xl0 20 cm~ 2 . 

The value determined by the ROSAT PSPC ob- 
servation is about 1.5(1.4-2.0) xl0 20 cm~ 2 (Nulsen and 
Bohringer 1995; Bohringer 1999). The observed HI col- 
umn density of our Galaxy by Stark et al. (1992) 
is 2.5xl0 20 cm~ 2 , which is slightly reduced to 1.8 x 



r (arcmin) 

Fig. 18. Projected line brightness profile of Fe-K of the 
EMOS (open circles) and EPN (crosses). Errors corre- 
spond to 68% confidence level. Solid line and dashed lines 
represent the results using APEC model and deprojected 
temperature profile, £;T w hoi e an( i ^Thardj respectively. We 
also plotted the profile for temperatures shifted by 30% 
(dotted line) 
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Fig. 19. The radial profile of hydrogen column density 
of EPN (crosses) and EMOS (diamonds) using the sin- 
gle temperature (gray) and the two temperature (black) 
MEKAL model fit. Holizontal line corresponds to the col- 
umn density obtained from Nulsen & Bohringer (1995). 



10 20 cm -2 , by Lieu et al. (1996). The observed A^h is con- 
sistent with the value obtained by ROSAT and also the 
Galactic column density. In addition, the power-law com- 
ponent from the central AGN shows no significant excess 
absorption. Thus, the intrinsic absorption should be neg- 
ligible in the M 87 system. 
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Table 6. Electron density profile 
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Table 7. result of double /3 fit 



7. Density profile 

We have assumed Nh to be the average value of those ob- 
tained at R > 2', since the higher Njj at R < 2' from the 
single MEKAL model may be artificial. We have fixed the 
C, N, O abundance to 0.51(1 + (i?/0.76') 2 ) -0 09 , which is 
the best fit /3-model of the abundance profile of these el- 
ements. Since the cooler component is distributed around 
the radio structures, and the hotter component is almost 
spherically symmetric, the latter should dominate. Thus, 
within 0.35 to 1.4', we fitted the EPN spectra with a two 
component MEKAL model and adopted the electron den- 
sity of the hotter component. The electron density (n e ), 
profile derived is shown in Table 6 and Figure This re- 
sults are mostly consistent with those obtained by Nulscn 
and Bohringer (1995), when considering the difference in 
the assumption of the abundance. 

We fitted the profile with a sum of two models, a 
compact and an extended one, both centered on M 87. We 
fixed the ft of the extended component to be 0.47. which is 
the value obtained by ROSAT All Sky Survey (Bohringer 
et al. 1994). The obtained best-fit parameters are shown in 
Table 7. The two [3 components have distinct parameters, 
and cross over at R ~ 2', or ~ 10 kpc. 



8. Dark matter profile 

The total mass M(R) within a 3-dimensional radius R, 
assuming hydrostatic equilibrium and spherical symmetry, 
is given by (e.g Sarazin 



M{R) 



kTR /dlnn d\nT 



Gfirrip \d\nR d In R 



where m p is the proton mass, k is the Boltzmann constant, 
G is the constant of gravity, and /i ~ 0.63 is the mean 
particle mass in units of m p . 
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Fig. 20. Electron density profile. The solid line is the 
best-fit doublc-/3 model. 



Based on the deprojected temperature and density 
profile obtained in §3.2 and §5, we have calculated the 
gravitational mass of the M 87 system, together with the 
gas mass profile using the power-law relation in Figure 
3, and the double-/? relation approximated as in Table 7. 
The upper and lower limit of the mass profile is calculated 
considering the uncertainty of the density gradient and 
the temperature. The temperature gradient is assumed to 
be smooth and expressed as the power-law relation. The 
upper and lower limit of the density gradient of the i-ih 
shell is obtained from the ratio of the value within i + 1-th 
shell to that within i — 1-th shell. As plotted in Figure pl| 
our result is mostly consistent with the dynamical mass 
profiles obtained by stars (Kronawitter et al. 2000) and 
globular clusters (Cohen & Ryzhov 1997). 

The calculated mass-to-blue light ratio (hereafter 
M/Lb) of the stellar population ranges between 6 and 10 
(e.g. Kronawitter et al. 2000). Assuming a constant stellar 
M/Lb of 8, in Figure |2l[ we also plotted the stellar mass 
profile, calculated from the luminosity profile in Giraud 
(1999). Within nearly 0.5r e , the stellar velocity dispersion 
profile indicates that the dynamical M/Lb is ~H (e.g. 
Kronawitter et al. 2000). Here, r e is the effective radius of 
M 87, which is 7.8 kpc. Our mass profile indicates that at 
R=l-2r e , the M/Lb is about 17 and beyond this radius, 
the slope of the mass profile increases and M/Lb reaches 
~ 200 at 10r e . 



9. Application of the cooling flow model 

9.1. Mass flow profile from imaging data 

The steady state cooling flow structure in a spherically 
symmetric cluster is obtained by, 
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Fig. 21. (upper panel) Integrated mass profiles of the to- 
tal gravitating matter (dashed lines for upper and lower 
limit and solid line for the double f3 model fit), the X- 
ray emitting plasma ( dot-dashed line), and the stellar 
component (dotted line) assuming stellar M / L-q to be 
8, in the M87 system. The dynamical mass obtained by 
stars (cross; Kronawitter et al. 2000) and globular clus- 
ters (open circles; Cohen & Ryzhov 1997) are also shown, 
(lower panel) Total mass-to-light ratio (M/L-q) of the sys- 
tem. The meanings of symbols and lines are the same as 
in the upper panel. 



M = 4irR 2 p g v 

(e.g. Sarazin 1988), where p g , v, and P are, respectively, 
the gas density, velocity, and pressure, and <f)(R) is the 
gravitational potential, A is the emissivity. 

Even within the optical radii of elliptical galaxies, the 
X-ray luminosity is a factor of 20 larger than heating by 
stellar mass loss (Matsushita 2001). Therefore, we can ne- 
glect this component as a sufficient heat source to balance 
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Fig. 22. Cooling time (upper panel) and M\ profile (lower 
panel). Data corresponds to Ms, assuming cooling to 
keV. The dotted line is the differential of the Mi, the mass 
deposition rate in each shell. 

cooling. The cooling time is lower than the Hubble time 
within R ~ 50 kpc (Figure |||). As shown in Fig. |2^, the 
Mi obtained is roughly proportional to the radius. At R ~ 
50 kpc, Mi is 10 M , which is consistent to the classical 
value obtained by Stewart et al. (1984) from the Einstein 
observation. 

9.2. Spectral fitting using a spectral cooling flow model 

With the RGS we can observe only the very central core of 
cooling flow clusters. In contrast, using the EPN, we can 
accumulate spectra within the entire cooling flow region. 

At R = 4', the derived cooling time is only a few Gyr 
and Mi is 8 M Q y _1 , respectively. The net effect of reso- 
nance scattering is negligible for the radiative output in 
this integrated circular region. Without any heating, the 
gas should cool down to keV, emitting X-rays at in- 
termediate temperatures. Therefore, the spectral cooling 
rate, Mg, within the radius should be 8AfQy _1 . 
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We fitted the deprojected spectra within R < 4' with 
a cooling flow model using the MEKAL code modified by 
photoelectric absorption. As in other cooling flow clusters 
observed by RGS, in order to fit the spectra, we still need 
a cut-off temperature, lowT, which is found to be 1.4 keV 
(Table 8). Figure^ shows the EPN spectrum within R < 
4', fitted with the cooling flow model. We also plotted 
in the figure, another cooling flow model with lowT=0.1 
keV, fitted only to the data in the spectral range above 
1.2 keV. The difference in the continuum level between 
0.2 and 0.5 keV is due to the instrumental low-energy tail 
of Fe-L and O-K lines. The component cooling radiatively 
to keV should emit strong Fe-L lines between 0.6 to 1.0 
keV, which are not seen in the spectrum. 

When adding an isothermal MEKAL component, we 
can fit the spectra with a cooling flow component with 
lowT=0.1 keV. However, M s is only 0.8 M y _1 , which is 
a factor of 10 smaller than M\ at the same radius. 

9.2.1. Test of a bimodal abundance model 

Fabian et al. (2001) made a proposal to solve the cool- 
ing flow problem by a scenario in which metals in the 
ICM are not uniformly distributed. In this suggestion, a 
metal poor part of the gas cools without emitting lines 
and a metal rich part cools rapidly. In this way, the total 
strength of Fe-L line emission should be reduced. A bi- 
modal abundance model, with 0.01 solar abundances and 
10 times higher abundances than the single abundance 
model and mass deposition rates of 7.3 and 0.7 M & /yr, 
respectively, is shown in Figure 23. It is very different from 



the observed spectrum in the Fe-L structure and strength 
of K-q lines (Figure |23|), since it reduces not only Fe-L 
lines but also K lines of Si, S, Ar, Ca, and Fe (see also 
Bohringer et al. 2002). In order to obtain the observed 
strength of K lines with these mass deposition rates, the 
abundances in the high metallicity component should be 
~ 100 solar. Therefore, we have tried another bimodal 
metal abundance model, where the mass deposition rates 
of the two components are allowed to be free parameters, 
and fitted the spectra above 1 keV. Figure ^3] also shows 
the best fit model, extrapolated to lower energies. With 
the mass deposition rates of 4 and 4 M & /yr, for the metal 
rich and metal poor components respectively, the spectra 
above ~ 1 keV can be fitted. Because of the dependence 
of the cooling function on metallicity, the Fe-L strength 
below 0.9 keV is slightly reduced compared to that of the 
single abundance model. However, the model shows very 
strong Fc-L lines below 0.9 keV. We then added absorp- 
tion for the metal rich components and then allowed the 
abundances of the metal rich component to vary. Although 
including a column density of 5xl0 21 cm -2 dramatically 
improved the fit, the fit is much worse compared to the fit 
with the cut-off temperature. In addition, such high ab- 
sorption heavily reduces the spectra around the 0.67 keV 
O line. Increasing the O abundance of the metal rich com- 
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Fig. 23. Deprojected spectra within R < 4' fitted with a 
cooling flow model with a cut-off temperature, lowT=1.44 
keV (solid line, a). The bold dotted line (b) shows the 
model with cooling to 0.1 keV. The thin dashed lines of 
the upper panel, (cl) and (c2), corresponds to the model 
with a model with 10 times higher metallicity and a model 
with low metal abundance, respectively. The spectral cool- 
ing flow rates are summarized in the panel. The sum of 
(cl) and (c2) is also plotted. The dotted lines of the lower 
panel, (dl) and (d2) are components of another bimodal 
metal abundance mdoel, and the solid lines of (el) and 
(e2) corresponds to the model with excess absorption ap- 
plied for the metal rich components. 



ponents does not solve the problem, since it reduces the 
strength of Fe-L lines due to the difference in the cooling 
function. Therefore, there must exist an unabsorbed com- 
ponent whose O abundance is not so peculiar. For intrinsic 
absorption, more details are discussed also on Bohringer 
et al. (2002). 

In summary, any multi-abundance model cannot ex- 
plain the observed Fe-L profiles. In order to explain it, we 
need a sharp cut off in the temperature distribution. 
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Table 8. Result of spectral fitting using a spectral cooling flow model 



9.2.2. Radial profile of spectral cooling flow component 

The observed radial profile of Mi is not constant, but 
roughly proportional to R. The differential of Mi between 
a shell and next shell implies that we should observe a 
cooling flow component in the spectrum of each radial 
shell. We thus fitted deprojected spectra with a cool- 
ing flow model with lowT=0.1 keV, and an isothermal 
MEKAL model. The result is summarized in Table 8 and 
Figure ^2|. Within R < 1', the fit has improved from the 
single temperature fit. However, the Ms values, are a fac- 
tor of 10 to 20 lower than the differential of M\. Only for 
the innermost shell, Ms and M\ are consistent with each 
other. 

10. Discussion 

10.1. Temperature structure 

The deprojected spectra of M 87 strongly indicate that ex- 
cluding the region associated to the radio structures, the 
ICM is single phase; at a given radius, the ICM is domi- 
nated by one temperature component. We cannot detect 
any multi-phase cooling flow component as seen by Allen 
et al. (2001). The ICM is also not characterized by a two- 
phase structure as suggested by ASCA observations of the 
Centaurus cluster (Fukazawa et al. 1994; Ikebe et al. 1999), 
and M 87 (Matsumoto et al. 1996). 

When excluding the regions with the radio struc- 
tures, the ICM is nearly spherically symmetric. Using the 
MEKAL model, a single phase model can fit the depro- 
jected spectra well at R > 2'. The temperature profile 
has a positive temperature gradient, ~ 1.7 keV at R = V 
and ~ 2.5 keV at R=15'. The consistency of the temper- 
atures obtained by the continuum, Fe-L region, Si, and 
S line ratios, Fe-L to Fe-K ratio indicate that the ICM 
is dominated by a single temperature component. The re- 
gions associated to the radio structures have an additional 
temperature components of 1 keV. At R < 2', since it is 
difficult to filter the complicated radio structures within 
2' with the XMM spatial resolution, we observe the major 



ICM components with the temperature of 1.7 keV, and 
small amount of the 1 keV components. The spatial dis- 
tribution of the latter component indicates that the cooler 
component also relates to the region corresponding to the 
radio lobes (see also Belsole et al. 2001). 

The R-S model gives a temperature lower by 20% than 
MEKAL model. In order to fit the spectra, the R-S model 
requires an additional high temperature component. How- 
ever, the two temperature R-S model cannot explain the 
observed Si and S line ratios. 

Except for the Fe-K line strength, the APEC temper- 
atures are consistent with the MEKAL models within 10- 
15%. Due to changes of the strength of Kq lines, the single 
phase APEC model cannot explain the strength of Fe-K 
lines, although it can fit the Fe-L and continuum simulta- 
neously. The S line ratio is discrepancy by a factor of 1.5, 
but it introduces only a 10% difference in the temperature 
structure. 

Resonant scattering of line emission should be impor- 
tant in dense core of galaxy clusters. (Gil'fanov et al. 1987, 
Tawara et al. 1997). Shigeyama (1998) calculated the ef- 
fect on M 87 system, and found that surface brightness 
profile considerably decreased within 1-2' because many 
Fe-L lines are optically thick. In M87, Bohringer et al. 
(2001) suggested that the resonant scattering may explain 
the observed sharp abundance drop at center. This effect 
may also change the spectra and the Fe-L line profile in the 
central region. Part of the 1 keV component may be due 
to this effect. Therefore, in order to study the temperature 
structure within 1-2' arcmin exactly, we have to consider 
effect of the resonance scattering. However, since the tem- 
perature obtained from the continuum spectra is close to 
that from the whole energy band, the main temperature 
component should be similar to the obtained values. At 
least using the MEKAL model, the strength of the Fe-L 
and Fe-K are consistent with each other, which have a 
different dependence on the resonance line scattering. In 
addition, most of the 1 keV component should be related 
to the radio activity which has a complicated spatial distri- 
bution (Belsole et al. 2001). Most of the resonantly scat- 
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tered light will be reemitted at slightly larger projected 
radii. Therefore integrating over the region, R < 4' will 
reduce the net effect of scattering and the spectra will 
hardly be disturbed. 

Matsushita (2001) found that normal elliptical galax- 
ies have /3 S pec=l- The central stellar velocity dispersion of 
M 87 is 350 km/s (Whitemore et al. 1985). The assump- 
tion of /3 sp cc=l implies a temperature of 0.8 keV, which 
agrees well with the observed central ICM temperature of 
M 87. The ICM temperature at 10-15' of 2.5 keV also 
agrees with the velocity dispersion of early-type galaxies 
of the Virgo cluster, which is 570 km/s (Binggeli and Tam- 
mann 1987). The observed constancy in the temperature 
suggests that the central cool temperature mainly reflects 
the galaxy potential, even though some part of the cool 
component may be produced by a small cooling flow in 
the central region. 

10.2. Intrinsic absorption 

The hydrogen column density is almost flat within the field 
of view. In addition, the spectrum of the central AGN does 
not show any absorption feature. Therefore, intrinsic ab- 
sorption should be negligible in the M 87 system. If there 
is excess Njj of the order of 10 20 cm~ 2 and the absorbing 
matter is uniformly distributed within 10 kpc and 50 kpc, 
the mass of the absorbing matters is 10 8 M Q and 10 9 M Q , 
respectively, which is much smaller than the cooling flow 
mass deposition rate accumulated over a Hubble time. 

10.3. Potential structure of ellipticals and clusters of 
galaxies 

We have determined the gravitational mass profile of the 
M 87 system, directly from the observed density and tem- 
perature profiles. At 0.5-2 r e , M/Lb of the system is con- 
stant at ~ 17. Beyond this radius, it starts to increase to 
- 200 at 80 kpc. 

This structure is essentially the same feature as seen in 
the Fornax cluster (Ikebe et al. 1996), the Centaurus clus- 
ter (Ikebe et al. 1999) and NGC 4636 (Matsushita et al. 
1998). Figure [24| summarizes the gravitational mass pro- 
files of the 3 cD galaxies and NGC 4636. These galaxies 
have similar hierarchical mass profiles; the gravitational 
mass gradually increases up to several r e , where the grav- 
itational to stellar mass ratio is common among galaxies. 
Then, M/L start to increase steeply. 

Before the total to stellar mass ratio increases steeply, 
the total M/Lb is ~ 17 at 1 r e , and ~ 30 even at 10 
r e . This feature should reflect the potential of elliptical 
galaxies. The M/Lb of the central core of these galaxies 
obtained by the dynamics and the same ratio obtained by 
the stellar population synthesis models are consistent with 
each other giving a typical value of 6-10 (e.g. Kronawitter 
et al. 2000). Considering the observed metallicity gradi- 
ent of stars, the stellar M/Lb should be slightly smaller 
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Fig. 24. Integrated mass profiles of the total gravitat- 
ing matter scaled by Lb, and the stellar component (solid 
line) in the M 87 (solid line; L x / L a =22), the Centaurus 
cluster (dotted line; Ikebe et al. 1999; L x /L a =A7), NGC 
1399 (dashed line; Ikebe et al. 1996 ;L X /L (T =3), and NGC 
4636 (gray line; Matsushita et al. 1998; L x /L a =Yl). L x 
and L a means the ISM luminosity within 4r e and kine- 
matical energy input from stellar mass loss (Matsushita 
2001). 

in the outer regions (e.g. Kodama & Arimoto 1997). Our 
result suggests that the dark mass distribution of elliptical 
galaxies is slightly more extended than those of stars. It 
is thought that elliptical galaxies have a 10 times larger 
amount of dark mass than stellar mass (e.g. Ciotti et al. 
1991). However, the galaxies themselves may not contain 
such large amount of mass. This result should be impor- 
tant for the study of the origin of the dark matter content 
in elliptical galaxies and also for the study of the formation 
and evolution of these galaxies. 

The scatter of the X-ray to blue luminosity ratio 
of early-type galaxies has long been a problem (e.g. 
Canizares et al. 1987). Matsushita (2001) found that L x 
of most of the early-type galaxies is well explained by a 
kinematical heating of the gas supplied by stellar mass 
loss {L a )- All of the galaxies which have much larger L x 
than L a have very extended X-ray emission. Figure [24| 
shows that the galaxies with a larger L x / L a have smaller 
radii of the starting point of the strong increase of M/L, 
when scaled r e . In other words, the X-ray luminosity of cD 
galaxies may be determined in relation to their potential 
structure, as indicated by Matsushita (2001). 

10.4- Problem of standard cooling flow model 

The standard cooling flow model indicates that the ICM is 
multi-temperature on small scales, since it is thought that 
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mass is deposited within the whole cooling flow region. 
The cooling matter should emit strong Fe-L lines below 
0.9 keV. We cannot detect a significant spectral cooling 
flow component from the whole field of view of the detec- 
tor and the ICM at any given radius is dominated by a 
single temperature component. The upper limit on Ms is 
an order of magnitude smaller than Mi . A bimodal metal 
abundance distribution model cannot explain the observed 
Fe-L profile. 

The resonance scattering may reduce the strength of 
some resonance lines. However, all the Fe-L lines below 
0.9 keV, which includes both resonance lines and non- 
resonance lines, are suppressed compared with the ex- 
pected value for a cooling flow. In addition, the spectral 
cooling flow component is also small where the resonance 
scattering is not effective. Therefore, the resonant scatter- 
ing can not be responsible for the suppresion of the Fe-L 
lines below 0.9 keV. 

At 1 r e (1.6') of M 87, M\ is ~ 4M yr" 1 . In contrast, 
the upper limit on Ms is ~ 0.4M Q ?/r~ :L . This value is 
close to the stellar mass loss rate within 1 r e , which is 0.5 
Moyr" 1 . At this radius, i coo i is 10 9 yr and the gas mass is 
2 10 9 M Q . Therefore, considering the observed abundance 
gradient, the standard cooling flow model indicate that 
the gas should be diluted by metal poor surrounding gas. 
However, the observed Fe abundance of M 87 within lr e , 
~ 1 solar, is similar to those of normal X-ray luminous 
galaxies, which are also about ~ 1 solar (Matsushita et 
al. 1997; 2000). Note that the observed ISM abundance of 
NGC 4636 from the RGS spectra is 0.9 solar (Xu et al. 
2001). In these systems, the enrichment has to be com- 
pared with a mass flow rate of ~ lM Q yr _1 as implied by 
the standard cooling flow models for these galaxies. For 
the case of NGC 4636, the mass flow rate is determined 
to be ~ lMoyr" 1 (Bregman, Miller, & Irwin 2001). The 
effect of the dilution by ICM should be determined by 
a ratio of the mass flow rate and the optical luminosity 
since metals ejected from a galaxy should be proportional 
to the stellar luminosity. Within lr e , the former rate of 
NGC 4636 is a factor of 10 smaller than that of M 87, 
while the latter is only a factor of 2.5 smaller. In a much 
larger cooling flow as exists in M 87 in which the Fe en- 
richment is supplied by a stellar population comparable 
to that of normal ellipticals, the relative enrichment will 
be less and the abundance gradient must be reduced in 
contrast to the observed data. Therefore, the overall cool- 
ing flow expected from the classical galaxy cluster cooling 
flow model, should not exist in M 87, although a small 
scale cooling flow, with mass deposition rate close to the 
stellar mass loss rate, may exists in the central core of the 
classical cooling flow regions. 

11. Conclusion 

The major result of this paper is strong evidence that the 
ICM in the halo of M 87 is isothermal locally, except prob- 



ably for the very central region, where more than one tem- 
perature component is present. This could be partly due 
to the interaction effects of the ICM with the jet and the 
radio lobes of M 87. This finding and the fact that no spec- 
tral signature of low temperature components below tem- 
peratures of 0.8 keV are observed is in disagreement with 
the standard cooling flow model which predicts a strongly 
multi-phase structure of the ICM and a distinct tempera- 
ture distribution and resulting spectrum with clear signa- 
ture of low temperature components. We have also shown 
that strong inhomogeneities in the metal distribution can- 
not resolve this problem. Therefore we conclude that the 
scenario for the dense gas with short cooling time in the 
centers of clusters needs to be revised. Probably heating- 
processes that can substantially reduced the mass deposi- 
tion have to be reconsidered, like the energy input of the 
central AGN. The central temperature closely reflects the 
gravitational potential depth of the central galaxy, rather 
than the existence of a cooling flow. 
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